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This  study  investigates  the  dynamic  relationship  between  coal  consumption,  economic  growth,  trade 
openness  and  C02  emissions  in  case  of  India.  In  doing  so,  Narayan  and  Popp,  Journal  of  Applied  Statistics 
2010;  37:1425-1438,  structural  break  unit  test  is  applied  to  test  the  order  of  integration  of  the  variables. 
Long  run  relationship  between  the  variables  is  tested  by  applying  the  ARDL  bounds  testing  approach  to 
cointegration  developed  by  Pesaran  et  al.  Journal  of  Applied  Econometrics  2001;  16:289-326. 

The  results  confirm  the  existence  of  cointegration  for  long  run  between  coal  consumption,  economic 
growth,  trade  openness  and  C02  emissions.  Our  empirical  exercise  indicates  the  presence  of  environ¬ 
mental  Kuznets  curve  (EKC)  in  long  run  as  well  as  in  short  run.  Coal  consumption  as  well  as  trade 
openness  contributes  to  C02  emissions.  The  causality  analysis  reports  the  feedback  hypothesis  between 
economic  growth  and  C02  emissions  and  same  inference  is  drawn  between  coal  consumption  and  C02 
emissions.  Moreover,  trade  openness  Granger  causes  economic  growth,  coal  consumption  and  C02 
emissions. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  Kyoto  Protocol  requires  that  industrialized  countries 
reduce  their  collective  emissions  of  greenhouse  gasses  by  5.2% 
of  1990  levels  by  the  period  2008-2012.  The  country-specific 
targets  in  the  Kyoto  Protocol  may  be  difficult  for  some  nations  to 
achieve.  Developing  countries,  including  India,  have  absolved  of 
any  responsibility  towards  reducing  emissions  in  the  first  com¬ 
mitment  period,  that  is,  2008-2012,  of  the  Kyoto  Protocol. 
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However,  India  estimates  a  24%  cut  by  2020  in  its  carbon  intensity 
i.e„  the  amount  of  carbon  dioxide  emitted  for  each  unit  of  GDP, 
compared  with  2005  levels,  and  by  2030,  it  estimates  it  could 
achieve  a  reduction  in  its  carbon  emissions  by  37%  from  2005 
levels,  according  to  provisional  government  figures  (The  Indian 
Express  [71]). 

Theoretically,  the  environmental  Kuznets  curve  (EKC)  hypothesis 
postulates  the  existence  of  an  inverted-U  shape  relationship 
between  real  GDP  per  capita  and  measures  of  environmental 
degradation  such  as  S02  and/or  C02  emissions.  However,  the 
empirical  evidences,  either  using  time  series  and/or  pooled  data  of 
a  group  of  countries,  on  the  EKC  hypothesis  vary  from  country  to 
countiy,  instead.  Further,  the  results  are  not  uniform  across  pollu¬ 
tants.  This  has  created  tow  problems  being  faced  by  environmental 
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policy  makers:  to  ensure  that  useful  knowledge  informs  policy 
(without  being  misused  and/or  distorted)  and  to  understand  how 
to  respond  to  this  knowledge  (Boehmer-Christiansen  [1]).  However, 
in  the  present  context  we  have  limited  ourselves  to  provide 
evidence  of  the  EKC  hypothesis  for  India.  Our  contribution  lies, 
particularly  in  two  directions.  First,  we  employed  the  ARDL  method 
which  is  amenable  for  short  time  series  data  as  in  this  paper  and 
second,  we  provide  empirical  evidence  of  the  EKC  by  incorporating 
coal  consumption  and  trade  as  additional  determinant  of  C02 
emissions  in  case  of  India. 

The  rational  for  selecting  India  for  our  analysis  is  that  it  has 
implemented  a  variety  of  programs  and  policy  initiatives  since 
the  introduction  of  the  National  Forest  Policy  in  1988.  India  has 
become  one  of  the  fast  growing  countries  next  to  China.  However, 
in  the  same  time  consumption  of  coal  has  increased  rapidly  which 
was  35.55  mote  in  1965  has  reached  the  height  of  249.86  mote  in 
2009.  The  close  relationship  between  coal  consumption  and 
economic  growth  of  India  is  evident  from  following  Fig.  1.  Further 
if  we  see  the  relationship  between  percentage  growth  rates  of 
GDP  and  coal  consumption,  we  find  that  recently  percentage 
growth  rate  in  the  coal  consumption  is  more  than  twice  of  the 
percentage  GDP  growth  rate.  This  situation  is  evident  in  Fig.  2. 

With  this  background  we  set  objective  to  test  the  environ¬ 
mental  Kuznets  curve  hypothesis  (EKC)  in  augmented  equation 
framework.  According  to  the  EKC  hypothesis  environmental 
degradation  increases  at  initial  level  of  economic  growth  and 
then  starts  to  decrease  at  a  higher  level  of  economic  growth. 
Hence,  the  relationship  between  measure  of  environmental 
degradation,  (in  our  case  it  is  C02  emissions)  and  measures  of 
economic  growth  (in  our  case  it  is  measured  by  real  GDP  per 
capita)  is  the  inverted-U  shaped  curve.  The  objective  of  present 
study  is  to  investigate  the  EKC  for  the  Indian  economy  over  the 


Fig.  1.  Coal  consumption  and  GDP  per  capita. 


period  of  1966-2011.  In  addition,  we  also  include  coal  consump¬ 
tion  and  trade  openness. 

The  rest  of  the  paper  consists  of:  Section  2  reviews  a  selected 
literature  encompassing  the  EKC  and  the  variables  listed  above, 
Section  3  has  the  theoretical  and  the  econometric  model  includ¬ 
ing  the  ARDL  estimation  strategy.  The  empirical  results  are 
reported  in  Section  4  followed  by  the  conclusion  and  policy 
implications. 


2.  Literature  review 

There  has  been  increasing  attention  on  the  impact  of  economic 
growth  on  environment  degradation  since  the  last  few  decades  of 
the  previous  century.  The  pioneering  attempt  in  this  area  was 
made  in  the  early  1990s  by  Grossman  and  Krueger  [2,3]  which 
investigated  the  environmental  impacts  of  the  North  American 
Free  Trade  Agreement.  Their  study  postulated,  estimated  and 
ascertained  an  inverted-U  shaped  relationship  between  measures 
of  several  pollutants  and  real  GDP  per  capita  which  was,  con¬ 
temporaneously,  confirmed  by  Shafik  and  Bandyopadhyay  [4]  and 
Panayotou  [5j.  Theoretically,  the  EKC  does  not  only  depend  on 
levels  of  per  capita  GDP,  but  also  on  a  series  of  factors  causing 
changes  in  economic  growth  that  can  affect  environment.  In 
general,  economists  analyze  mechanisms  behind  the  EKC  by 
examining  scale  effect,  structural  effect  and  technique  effect 
(Song  et  al.  [6]).  However,  few  scholars  (for  example,  Panayotou 
[5])  believe  that  EKC  is  caused  by  up-gradation  from  the  adjust¬ 
ment  of  economic  structure.  According  to  the  structural  effect 
hypothesis  (Stern  [7])  economic  development  passes  through  the 
various  states  starting  form  preliminary,  rapid-development  and 
high-grade,  industrial  structure  first  upgrades  from  agriculture  to 
a  high-pollution  industry,  finally  turns  to  information  concen¬ 
trated  industry,  which  leads  to  the  improvements  in  environ¬ 
mental  quality.  However,  as  Stokey  [8]  pointed  out,  due  to  the 
technique  effect  economic  growth  can  break  through  one  thresh¬ 
old  point  after  arriving  at  a  certain  stage  of  economic  develop¬ 
ment.  Hence,  at  a  low  income  level,  only  the  high  pollution 
technique  can  be  used,  but  once  leaping  over  the  threshold  point 
of  economic  development,  cleaner  technologies  can  be  adopted 
which  lowers  the  degradation  in  the  environmental  quality. 
Further,  some  scholars  attribute  the  demand  factors  to  the  cause 
of  EKC  (for  example,  Lopez  [9]),  which  asserts  that  demand  for  a 
clean  environment  will  be  increased  over  the  real  income  per 
capita.  Giving  the  importance  of  scale  effect,  Andreoni  and 
Levinson  [10]  suggested  that  in  the  static  model  of  single  depart¬ 
ment,  the  EKC  can  be  derived  technically,  only  if  pollution  control 
is  increasing  in  scale. 

Further,  Suri  and  Chapman  [11]  bring  the  contribution  of 
industrial  products  of  the  imports  and  exports  to  industrial 
products  of  national  production  into  the  analytical  framework 
of  the  EKC.  That  is  to  say,  the  low  emissions  correspond  to  the 
growth  of  industrial  products  of  the  imports,  while  the  high 
emissions  correspond  to  that  of  the  exports.  This  implies  that 
there  is  a  strong  relationship  between  trade  and  environmental 
quality,  so  the  evolution  of  environmental  quality  can  be  pre¬ 
dicted  effectively.  The  empirical  studies  of  EKC  started  by  Gross- 
man  and  Krueger  [2]  and  followed  by  Lucas  et  al.  [12],  Wyckoff 
and  Roop  [13],  Suri  and  Chapman  [11],  Heil  and  Selden  [14],  Friedl 
and  Getzner  [15],  Stern  [7],  Nohman  and  Antrobus  [16],  Dinda 
and  Coondoo  [17]  and  Coondoo  and  Dinda  [18]  but  presented 
mixed  empirical  evidence  on  the  validity  of  EKC.  Song  et  al.  [6], 
Dhakal  [19],  Jalil  and  Mahmud  [20]  and,  Zhang  and  Cheng  [21] 
supported  the  existence  of  EKC  in  China.  The  findings  of  Fodha 
and  Zaghdoud  [22]  revealed  the  existence  of  EKC  between  the  S02 
emissions  and  economic  growth  but  not  for  the  C02  emissions  in 
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Tunisia.  In  contrast,  Akbostancl  et  al.  [23]  did  not  support  the 
existence  of  EKC  in  Turkey  and  same  conclusion  was  drawn  by 
Saboori  et  al.  [64]  for  Indonesia.  They  argued  that  C02  emissions 
are  automatically  reduced  due  to  the  rapid  pace  of  economic 
growth. 

Recent  literature  documented  alliance  of  economic  growth 
with  energy  consumption  and  environmental  pollution  to  inves¬ 
tigate  the  validity  of  EKC.  The  relationship  between  economic 
growth,  energy  consumption  and  C02  emissions  have  also  been 
researched  extensively  both  in  the  country  case  and  panel  studies. 
For  example,  Ang  [24]  found  stable  long  run  relationship  between 
economic  growth,  energy  consumption  and  C02  emissions  for 
French  economy  while  Ang  [25]  also  got  similar  results  for 
Malaysia.  Ang  [24]  showed  that  causality  is  running  from  eco¬ 
nomic  growth  to  energy  consumption  and  C02  emissions  in  the 
long  run  but  energy  consumption  Granger-causes  economic 
growth  in  the  short  run.  In  case  of  Malaysia,  Ang  [25]  reported 
that  output  increases  C02  emissions  and  energy  consumption. 
Ghosh  [26]  documented  that  no  long  run  causality  between 
economic  growth  and  C02  emissions  and  in  the  short  run, 
bidirectional  causality  exists  in  India.  On  contrary,  Tiwari  [66] 
noted  unidirectional  causality  running  from  C02  emissions  to 
economic  growth  and  C02  emissions  and  energy  consumption 
and  C02  emissions  Granger  cause  each  other.  In  case  of  Indonesia, 
Hwang  and  Yoo  [69]  noted  bidirectional  causality  between  energy 
consumption  and  C02  emissions  in  long  run  but  in  short  run, 
energy  consumption  and  C02  emissions  are  cause  of  economic 
growth.  Latter  on,  Uddin  et  al.  [70]  investigated  the  relationship 
between  energy  consumption,  economic  growth,  trade  openness 
and  C02  emissions  in  case  of  Srilanka.  Their  results  found  that 
economic  growth  Granger  causes  energy  consumption  and  C02 
emissions. 

For  the  panel  studies,  Apergis  and  Payne  [27]  investigated  this 
relationship  for  six  Central  American  economies  using  panel 
VECM.  It  is  evident  that  energy  consumption  is  positively  linked 
with  C02  emissions  and  EKC  hypothesis  has  been  confirmed.  Lean 
and  Smyth  [28]  and  Apergis  and  Payne  [29]  reached  the  same 
conclusion  for  the  case  of  ASEAN  countries  and  Commonwealth  of 
Independent  States,  respectively.  Narayan  and  Narayan’s  [30] 
empirical  evidence  also  validates  the  EKC  hypothesis  for  43  low 
income  countries.  In  addition,  Lean  and  Smyth  [28]  noted  a  long 
run  causality  running  from  energy  consumption  and  C02  emis¬ 
sions  to  economic  growth  but  in  short  span  of  time,  energy 
consumption  Granger  causes  C02  emissions.  On  the  other  hand, 
Apergis  and  Payne  [29]  found  that  energy  consumption  and 
economic  growth  Granger  causes  C02  emissions  while  bidirec¬ 
tional  causality  is  found  between  energy  consumption  and 
economic  growth;  and  between  energy  consumption  and  C02 
emissions.  Chen  [31  ]  explored  this  issue  to  Chinese  provinces  and 
documented  that  industrial  sector’s  development  is  linked  with 
increase  of  C02  emissions  due  to  energy  consumption.  Similarly, 
Shiyi  [32]  investigated  the  relationship  between  industrial  sec¬ 
tor’s  development  and  C02  emissions  using  Chinese  provincial 
data  and  concluded  that  industrial  development  increases  C02 
emissions.  Pao  and  Tsai  [33],  Ozturk  and  Acaravci  [34]  and, 
Acaravci  and  Ozturk  [35]  also  validated  the  existence  of  the  EKC 
in  case  of  BRIC,  Turkey  and  Demark  and  Italy,  respectively.  In  case 
of  Pakistan,  Nasir  and  Rehman  [36]  and,  Shahbaz  et  al.  [37] 
reported  that  the  EKC  exists  while  a  rise  in  economic  growth 
and  energy  consumption  is  linked  with  an  increase  in  energy 
pollutants.  Iwata  [38]  investigated  the  empirical  existence  of  the 
EKC  in  28  countries  adding  Pakistan.  They  corroborated  the 
existence  of  the  EKC  by  nuclear  energy  consumption  taking  into 
account.  Saboori  et  al.  [65]  unveiled  that  long  run  as  well  as  short 
run  EKC  exists  in  case  of  Malasia.  Unidirectional  causality  from 
economic  growth  to  C02  emissions  further  malidates  the 


existence  of  KEC  for  Malaysian  economy.  Yeh  [67]  used  the 
quintiles  regressions  using  the  data  of  developed  and  developed 
countries  and  reported  the  existence  of  EKC  by  controlling  other 
macroeconomic  variables.  In  case  of  Romania,  Shahbaz  et  al.  [68] 
confirmed  long  run  relationship  between  economic  growth, 
energy  consumption  and  energy  pollutants.  Their  empirical  evi¬ 
dence  also  found  that  environmental  Kuznets  curve  (EKC)  exist 
both  in  long-and-short  runs.  Moreover,  energy  consumption  is 
major  contributor  to  energy  pollutants.  Democratic  regime  shows 
her  significant  contribution  to  decline  C02  emissions  through 
effective  implementation  of  economic  policies  and  financial 
development  improves  environment  i.e„  reduces  C02  emissions 
by  redirecting  the  resources  to  environment  friendly  projects. 

The  relationship  between  trade  openness  and  environmental 
degradation  has  also  been  investigated  empirically.  Grossman  and 
Krueger  [2])  argued  that  environmental  effect  of  international 
trade  depends  on  the  policies  of  an  economy.  There  are  two 
schools  of  thought  about  the  impact  of  international  trade  on  C02 
emissions.  One  argued  that  trade  openness  provides  an  offer  to 
each  country  to  have  access  to  international  market  which 
enhances  the  market  share  among  countries.  This  leads  the 
competition  among  the  countries  and  increases  the  efficiency  of 
using  scarce  resources  and  encourages  importing  cleaner  tech¬ 
nologies  to  decline  C02  emissions  (e.g.,  Runge  [39]  and  Helpman 
[40]).  Other  group  probed  that  natural  resources  are  depleted  due 
to  international  trade.  This  depletion  of  natural  resources  raises 
C02  emissions  and  causes  environment  quality  worsened  (e.g., 
Schmalensee  et  al.  [41  ];  Copeland  and  Taylor  [42];  Chaudhuri  and 
Pfaff  [43]).  In  country  case  studies,  Machado  [44]  indicated 
positive  link  between  foreign  trade  and  C02  emissions  in  Brazil. 
Mongelli  et  al.  [45]  concluded  that  pollution  haven  hypothesis  is 
existed  in  Italy.  Halicioglu  [46]  added  trade  openness  to  explore 
the  relationship  between  economic  growth,  C02  emissions  and 
energy  consumption  for  Turkey.  The  result  showed  that  trade 
openness  is  one  of  main  contributor  to  economic  growth  while 
income  raises  the  levels  of  C02  emissions.  Finally,  Jayanthaku- 
maran  et  al.  [47]  probed  the  relationship  between  economic 
growth,  energy  consumption,  trade  openness  and  C02  emissions 
in  case  of  India  and  China.  Their  empirical  evidence  confirmed  the 
validation  of  EKC  in  both  countries. 

3.  Modeling,  empirical  strategy  and  data  collection 

The  dynamic  relationship  between  income,  coal  consumption 
and  C02  emissions  is  investigated  by  including  trade  openness  as 
potential  determinant  to  C02  emissions.  For  this  purpose,  we 
followed  the  empirical  model  applied  by  Ang  [25,26]  for  France 
and  Malaysia,  Soytas  et  al.  [48]  for  United  States,  Jalil  and 
Mahmud  [20]  for  China,  Halicioglu  [46]  for  Turkey,  and  Shahbaz 
et  al.  [37]  for  Pakistan.  We  use  coal  consumption  as  an  indicator 
of  energy  rather  than  energy  consumption.  More  than  50%  energy 
demand  in  India  is  fulfilled  from  coal  consumption.  Following 
Shahbaz  et  al.  [37],  we  converted  all  the  series  into  natural 
logarithms  to  obtain  efficient  and  consistent  results.  The  relation 
is  specified  as  follows: 

CO  2,t=f(Yt,YlCt,TRt)  (1) 

In  the  log-linear  form  the  specification  is  written  as: 

lnC02,t  =  +  02lnYt  +  </>3lnY2  +(/)4lnCc+4>5\nTRc+^i[  (2) 

where,  C02  denotes  carbon  emissions  per  capita  (in  kt);  C  is  coal 
consumption  per  capita;  Y  and  (Y2)  mentions  real  GDP  per  capita 
and  its  square,  respectively.  TR  is  trade  openness  [(exports  +  im¬ 
ports  )/GDP)[  per  capita  /(  is  a  random  error  term.  Following  the 
EKC  hypothesis  we  expect  that  ip2>  0  and  tp3  <  0.  The  use  of  coal 
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is  detrimental  for  environment  and  it  is  expected  that  </>4  >  0.  The 
impact  of  trade  openness  on  C02  emissions  can  be  positive  or 
negative.  This  implies  that  <fi5<  0  if  government  enforces  envir¬ 
onmental  laws,  and  makes  possibility  to  import  of  environment 
friendly  capital  and  technology  to  be  used  in  production  process. 
But  Grossman  and  Krueger  [2,3]  argued  that  if  emissions  might  be 
generated  due  to  relocation  of  polluting  industries  from  devel¬ 
oped  economies,  a  practice  known  in  the  literature  as  the  'safe- 
haven  hypotheses’  the  </>5  >  0. 

We  applied  a  recent  unit  root  test  developed  by  Narayan 
and  Popp  [49]  to  attain  efficient  and  consistent  empirical  evi¬ 
dence  about  the  integration  order  of  the  variables.  This  test 
incorporates  about  two  structural  breaks  stemming  in  the  series 
which  is  based  on  two  regression  equations  of  the  form  as 
follows: 

yr  =  m-1+a  1  +P*t+0iD(T'B)u+e2D(rB)2t 

+  <51D(J'l  t_|  +  d2DU2t_]  +  J2  PM-j  +  et  (3) 

j  =  i 

y“2  =  pyt-i P*t+KiD(T'B)  1i:+k2d(t'b)2j 

T <51  DL/,  t  j  +  d2DU2  +  +y2DT2t_  ] 

k 

+  Y2PjAyt- j+et  (4) 

j= i 

where  et  is  a  white  noise  process,  such  that  et  ~  NIID(0,a2), 
D(rB\t  =  -l(t  =  TBi  +  -l),  i  =  1,2,  DU'i[  =  l(t>T'Bi),  DT'it  = 
1  (t  >T'Bi)(t-DT'Bi),  T'Bi,  i=l,  2,  denote  the  true  break  dates.  Eqs. 
(3)  and  (4)  are  IO-type  test  regression  for  Ml  and  M2,  respec¬ 
tively.  The  unit  root  null  hypothesis  of  p  =  l  is  tested  against  the 
alternative  hypothesis  of  p  <  1  by  using  the  t-statistics  of  p.  Since 
it  is  assumed  that  true  break  dates  are  unknown,  VBi  in  Eqs. 
(3)  and  (4)  has  to  be  substituted  by  their  estimates  tBi,  i=  1,  2,  in 
order  to  conduct  the  unit  root  test.  We  used  a  sequential 
procedure  following  Narayan  and  Popp  [49]  to  attain  fBi,  i=l,  2. 

We  apply  the  autoregressive  distributed  lag  model  or  the  ARDL 
bounds  testing  approach  to  cointegration.  The  ARDL  bounds 
testing  approach  to  cointegration  is  preferred  due  to  its  certain 
advantages  over  traditional  cointegration  approaches.  For  exam¬ 
ple,  the  ARDL  bounds  testing  is  flexible  regarding  the  integrating 
order  of  the  variables  whether  variables  are  found  to  be  stationary 
at  1(1)  or  1(0)  or  I(1)/I(0)\  The  Monte  Carlo  investigation  shows 
that  this  approach  is  superior  and  provides  consistent  results  for 
small  sample  (Pesaran  and  Shin  [50]).  Moreover,  a  dynamic 
unrestricted  error  correction  model  (UECM)  can  be  derived  from 
the  ARDL  bounds  testing  through  a  simple  linear  transformation. 
The  UECM  integrates  the  short  run  dynamics  with  the  long  run 
equilibrium  without  losing  any  information  for  long  run.  The 
empirical  formulation  of  the  ARDL  bounds  testing  approach  to 
cointegration  is  given  below: 

AlnCOj  =  (Xi  +aiT+aco2lnCOj!_1 

+  aylnYt_i  -Ft^lnyjLj  -t-aclnC^  +aT)jln  TRt„i 

P  9  r 

+  E  (XjAlnCO,  ]  +  E  «jAlnYt_j  +  E  akA\nY^ 
i  =  1  ]  =  o  k  =  o 

S  t 

+  ^  ^  a/AlnCt_j  H-  ^  ^  ocm Ain  + /if  (5) 

1  =  0  m  =  0 

AlnYf  =  ot-i  +  octT + aC02  lnCOf_i  +  ay  lnYt_i 


1  The  ARDL  approach  to  cointegration  is  applicable  if  variables  are  integrated 
at  1(1)  or  1(0)  or  I(1)/I(0). 


P 

+  $y2lnYf_1  +aclnCf_i  TR[-\  +  J2  PiAinYt~i 

i  =  1 

+  J2  ftAlnCO?.,  +  J2  AAln Y2 *C_,  +  E  AAlnCt_, 

j  =  0  k  =  0  1  =  0 

t 

+  An  Ain  TRt_i  +  fit  (6) 

m  =  0 

AlnYf  =  ai  +a7r+aco2lnCOf_1  +aylnYt_i 

p 

+  ay2lnYf_1  +aclnCf_i  +ocTR\n  TRt_ ^  ^  ft  Ain  Y2t_A 

i  =  1 

+  ^2  ftAlnCOf_1  +  ^feAlnYt_fe  +  ftAlnCt_/ 
j  =  o  k  =  o  /  =  o 

t 

+  ^  ^  f3m  Ain  TR[_m  +  /if  (7) 

m  =  0 

AlnCf  =  ol\  +a7T+ac02lnCOf_1  +  aylnYt_i 

p 

+  ay2lnYf_1  +aclnCt_i  TRt_ i  +  ^  ft-AlnCf.,- 

i  =  1 

+  ftAlnCOt,  +  J2  &k AlnYf_k  +  £  ftAln Y2t_x 

j  =  0  k  =  0  1  =  0 

t 

+  y  )  i9m Ain  TRt-m“F //f  (8) 

m  =  0 

Ain  TRC  =  ar+otTT+ac02lnC02r  +  aylnYt_i  -pa^lnY^j 

p 

-FaclnC^,  +areln  TRt_i  +  E  P.Aln  TRr_,- 

i  =  1 

+  E  PjAlnCOf  ,  +  E  PkAlnYt-k+  E  P/AlnY(_i 

j=0  k= 0  (=0 

f 

+  E!  Pin AhlCt  nil  fit  (9) 

m  =  0 

here,  to  test  for  the  existence  of  cointegration  calculated  F- 
statistic  is  used  to  compare  with  critical  bounds  generated  by 
Pesaran  et  al.  [51].  Two  critical  bounds  namely  upper  critical 
bound  (UCB)  and  lower  critical  bound  (LCB)  have  been  developed 
by  Pesaran  et  al.  [51].  In  our  specification  of  Eqs.  (5)-(9),  the  null 
hypothesis  which  is  tested  i.e.,  H0  :  aco2  =aY  =  <xY2  =xc  =  &tr  =  0 
against  alternate  hypothesis  H0  :  aC02  ^  aY  aY2  ^  ac  ¥=  otTR  ^  0  to 
test  for  cointegration.  The  F-test  is  non-standard  and  we  may  use 
the  LCB  and  UCB  developed  by  Pesaran  et  al.  [51  p  to  test  for 
cointegration.  If  computed  F-statistic  is  more  than  UCB,  there  is 
cointegration  between  the  variables.  If  computed  F-statistic  does 
not  exceed  the  LCB,  there  is  no  cointegration  relationship  and  if 
computed  F-statistic  falls  between  lower  and  upper  critical 
bounds  then  decision  regarding  cointegration  between  the  vari¬ 
ables  remains  uncertain3  .  However,  the  critical  bounds  generated 
by  Pesaran  et  al.  [51  ]  may  not  be  appropriate  for  small  sample  like 
ours  case  with  45  observations  hence,  we  used  critical  bounds 
developed  by  Narayan  [52].  Further,  we  have  conducted  the 
stability  tests  to  examine  the  stability  of  ARDL  bounds  testing 
estimates  by  utilizing  CUSUM  and  CUSUMSQtest  of  (Browns  et  al. 
[53]). 


2  Pesaran  et  al.  [51  ]  have  computed  two  asymptotic  critical  values — one  when 
the  variables  are  assumed  to  be  1(0)  and  the  other  when  the  variables  are  assumed 
to  be  1(1). 

3  In  such  case,  error  correction  method  is  appropriate  method  to  investigate 

the  cointegration  (Banerjee  et  al.  [54]).  This  indicates  that  error  correction  term 

will  be  a  useful  way  of  establishing  cointegration  between  the  variables. 


A.K  Tiwari  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  18  (2013)  519-527 


523 


After  confirming  the  evidence  of  cointegration  and  cointegra¬ 
tion  is  stable  we  moved  ahead  for  testing  the  causal  relationship 
between  the  test  variables  in  the  framework  of  vector  error 
correction  method  (VECM).  The  VECM  in  five  variables  case  can 
be  written  as  follows: 

l  m 

Ain  COjr  =  a0i  +  55  an  Ain  CO,_1  +  55  a22Aln  Yt_j 
i  =  1  i=  1 

n  o 

+  55  ^33  Ain  y2t- 1+  Y,  a44AlnCf_r 

k  =  1  r  =  1 

P 

+  OC55 Ain  +  ^|£CTf_i (10) 

s  =  i 

/  m 

AlnY  =  /?01  +  55  $11  Ain  Yt-i  +  55  feAln  COf_i 

i  =  i  j  =  i 

+  55  ^t-1  +  55  /^44^1nCt-r+  55  /^55 AlnTRf  ! 

fe  =  1  r  —  1  s=l 

+  t/2ECTt_i +/i2,-  (11) 

/  m 

Ain  Yj  =  0Oi  +  55  011  Ain  Yf_i  +  55  022Aln  COjLi 
i  =  1  j  =  1 

+  55  033 55  044^^nCt-l  +  55  055^1n^t-l 

k=l  r  —  1  s=l 

f/  jTCT,  ]  +  /i},  (12) 

l  m  n 

AlnCf  =  0O1  +  55  0n Ain Ct_i  +  55  022 AlnCO^j  +  55  033Aln Yt 

i  =  1  j  =  1  k  =  1 

0  p 

T  ^  ^  044AlnYf_1  4-  N  )  055AlnTRt^i  -3ij^ECl\  \  3- fi^ 

r  =  1  s  —  1 

(13) 

I  m 

AlnTkt  =  <5oi  +  55  <SiiAlnrRt_,  +  55  ^AlnCO^, 

i  =  1  j  =  1 

+  55  ^33AlnYt_t+  55  <544^1nYf_1  +  55  ^ss^lnQ-i 
fe  =  1  r  =  1  s=1 

-|  1/4 r.CT,  ]  ^ 4 ;  (14) 

where  uit  are  white  noise  residual  terms.  The  estimates  of  HCTt_  1 
also  shows  the  speed  of  convergence  from  short  run  towards  long 
run  equilibrium  path  in  all  models  depending  upon  the  sign  of  the 
coefficient  of  ECTt_  1  term.  The  usefulness  of  the  VECM  lies  in 
distinguishing  between  short-and-long  runs  causal  relationships 
in  one  hand  and  detecting  causality  in  long-run,  short-run  and 
joint  i.e.,  short-and-long  runs  on  the  other  hand.  The  significance 
of  ECTt- 1  can  be  tested  using  t-statistic  which  confirms  the  long- 
run  causal  relationship  whereas  the  short  run  causality  is 
detected  by  the  joint  y2  statistical  significance  of  the  estimates 
of  first  difference  lagged  independent  variables.  In  the  final  step 
we  used  Wald-test  to  test  the  joint  significance  of  estimates  of  the 
lagged  terms  of  independent  variables  and  error  correction  term 
which  confirms  the  existence  of  short-and-long  run  causality 
relations  (Shahbaz  et  al.  [55])  and  known  as  measure  of  strong 
Granger-causality  (Oh  and  Lee  [56]). 

The  data  on  these  variables  has  been  collected  from  World 
Development  Indicators  (CD-ROM,  2011).  Further,  World  Devel¬ 
opment  Indicates  has  been  used  to  collect  data  on  real  GDP  per 
capita,  coal  consumption  per  capita,  real  trade  (real  imports  per 
capita+real  exports  per  capita)  per  capita  and  C02  emissions  per 
capita.  The  study  covers  the  period  of  1966-2009. 


Table  1 

Ng-Perron  unit  root  test. 


Level 

Variables 

MZa 

MZt 

MSB 

MPT 

lnC02t 

-9.8577 

-2.2095 

0.2241 

9.2910 

lnY, 

-8.9336 

-1.8632 

0.2085 

11.0896 

InY? 

1.2003 

0.6869 

0.5722 

85.0399 

InC, 

-4.1594 

-1.1484 

0.2761 

19.0857 

InTRf 

-4.8699 

-1.4368 

0.2950 

18.0242 

1st  difference 

lnC02t 

-24.7353* 

-3.5163 

0.1421 

3.6862 

lnYt 

-19.2482** 

-3.1020 

0.1611 

4.7358 

InYf 

-23.9443* 

-3.4600 

0.1445 

3.8058 

lnCt 

-21.5425** 

-3.2815 

0.1523 

4.2324 

InTRf 

-22.4313* 

-3.3472 

0.1492 

4.0729 

*  Indicate  the  significance  at  the  1%  level. 
**  Indicate  the  significance  at  the  5%  level. 


4.  Empirical  results  and  discussions 

We  have  applied  the  ARDL  bounds  testing  approach  to  coin¬ 
tegration  to  find  long  run  relationship  between  economic  growth, 
coal  consumption,  trade  openness  and  C02  emissions  in  case  of 
India.  Although  bounds  testing  approach  is  flexible  regarding  the 
integrating  order  of  the  variables.  We  can  apply  it,  if  the  variables 
are  integrated  at  1(0),  1(1)  or  I(0)/I(1).  We  cannot  compute  the 
ARDL  F-statistics  for  long  run  relationship  if  any  variable  is  found 
to  be  stationary  at  1(2).  To  ensure  that  none  of  the  variables  is 
stationary  at  1(2)  or  beyond  that  order  of  integration,  we  applied 
Ng-Perron  [57]  unit  root  test.  The  results  are  reported  in  Table  1. 
This  empirical  evidence  indicates  that  all  the  series  have  unit  root 
at  their  level  form  but  found  to  be  stationary  at  1st  difference. 
This  shows  that  all  the  variables  have  1(1)  order  of  integration. 
The  empirical  evidence  provided  by  Ng-Perron  may  be  biased 
because  it  does  not  have  information  about  structural  break 
stemming  in  the  series. 

To  solve  this  issue,  we  applied  unit  root  test  with  two 
structural  breaks  developed  by  Narayan  and  Popp  [49].  The 
inappropriate  chose  of  structural  breaks  may  reject  null  hypoth¬ 
esis  and  provide  inconclusive  results.  The  technique  used  in  N-P 
unit  root  test  allows  in  determining  two  structural  breaks  in  the 
level  and  trend  of  the  series  endogenously.  The  unit  root  test 
developed  by  Narayan  and  Popp  [49]  indicates  two  structural 
breaks  efficiently  as  compared  to  other  structural  break  unit  root 
tests.  Table  2  reported  the  results  by  Narayan  and  Popp  [49]  unit 
root  test.  Our  empirical  exercise  reveals  that  all  the  variables  are 
found  to  be  non-stationary  (except  lnTRt,  which  is  in  model 
2  turns  to  be  stationary)  at  level  but  integrated  at  1st  difference. 
This  confirms  that  all  the  variables  are  integrated  at  1(1)  and 
robust4. 

Our  empirical  evidence  confirms  that  all  the  variables  are 
stationary  at  1st  difference  form  that  intends  to  apply  the  ARDL 
bounds  approach  in  examining  the  long  run  relationship  between 
the  series.  It  is  prerequisite  to  choose  appropriate  lag  length  of  the 
variables  to  capture  the  dynamic  relationship  between  the  vari¬ 
ables.  The  statistics  of  F-test  are  very  much  elusive  with  the 
selection  of  lag  order  of  the  variables.  Inappropriate  lag  length 
selection  provides  biased  results  and  in  resulting,  decision  about 
cointegration  is  worthless  (Liitkepohl  [58]).  We  have  used  A1C  and 
SBC  criterion  to  choose  lag  order  but  our  decision  is  based 
minimum  value  of  AIC.  Liitkepohl  [58]  suggested  that  lag  order 
selected  by  AIC  provides  efficient  and  consistent  results.  The 


4  The  results  of  N-P  at  1st  difference  are  not  provided  here  due  lack  of  space 
and  are  available  from  authors  upon  request. 
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Table  2 

Narayan  and  Popp  [49]  Structural  Break  Unit  Root  Test. 


Variable 

Model  1(M1) 

Model  2(M2) 

T-statistic 

TB1 

TB2 

T-statistic 

TB1 

TB2 

lnC02t 

-4.0378 

1976 

2000 

-4.3381 

1976 

2000 

InYf 

0.7720 

1974 

1978 

0.4483 

1974 

1978 

InY7 

1.0668 

1974 

1978 

0.7999 

1974 

1978 

lnCt 

-1.1955 

1973 

1989 

-1.6512 

1977 

1989 

InTRt 

-3.0225 

1974 

1981 

—  5.3848  * 

1984 

1994 

Note:  Model-1  assumes  two  breaks  in  level  and  Model-2  assumes  two  breaks  in 
level  as  well  as  slope. 

**  Denotes  significance  at  5%  level. 


Table  3 

Results  of  ARDL  cointegration  test. 


Variable 

lnC02f 

InY, 

InY7 

lnCt 

lnTR, 

F-statistics 

8.9129** 

7.1597*** 

7.1578“ 

7.4063*** 

5.3858 

Critical  values# 

1%  level 

5%  level 

10%  level* 

Lower  bounds 

10.150 

7.135 

5.950 

Upper  bounds 

11.130 

7.980 

6.680 

Diagnostic  tests 
R2 

0.8325 

0.9994 

0.9995 

0.6249 

0.8018 

Adj-R2 

0.6802 

0.9990 

0.9991 

0.3796 

0.6037 

Durbin- Watson 

1.6951 

2.2267 

2.2155 

2.0707 

2.3233 

*  Critical  values  bounds  are  from  Narayan  [52]  with  unrestricted  intercept  and 
unrestricted  trend. 

*  Show  the  significance  at  1%  level. 

**  Show  the  significance  at  5 %  level. 

***  Show  the  significance  at  10%  level. 


Table  4 

Results  of  Johansen  cointegration  test. 


Hypothesis  Trace  statistic  Maximum  eigen  value 


R= 0 

105.8022* 

46.7978* 

R  <  1 

59.0044* 

32.3417** 

R<  2 

26.6626 

17.6688 

R<  3 

8.9938 

8.6642 

R  <  4 

0.3295 

0.3295 

*  Show  significant  at  1%  level. 
**  Show  significant  at  5%  level. 


reason  is  that  A1C  has  superior  properties  and  high  explaining 
power  as  compared  to  SBC.  Our  empirical  evidence  shows  that  we 
cannot  take  lag  more  than  2  because  data  sample  consists  of  45 
years  Table  3.  We  have  used  yearly  data  over  the  period  of  1966- 
2009  of  Indian  economy5. 

After  choosing  lag  order,  next  step  is  to  compare  our  computed 
F-statistics  with  lower  and  upper  critical  bounds  generated  by 
Narayan  [52].  The  results  show  that  our  calculated  F-statistics 
exceed  upper  critical  bounds  at  5%  and  10%  levels,  respectively, 
once,  C02  emissions  per  capita,  real  GDP  per  capita  (squared  of 
real  GDP  per  capita)  and  coal  consumption  per  capita  treated  as 
respond  variables.  The  results  are  reported  in  Table  4. 

This  implies  that  three  cointegrating  vectors  are  found  recom¬ 
mending  the  long  run  relationship  between  coal  consumption, 
economic  growth,  trade  openness  and  C02  emissions  in  case  of 
India  over  the  study  period  of  1966-2011. 

The  robustness  of  long  run  results  established  by  the  ARDL 
bounds  testing  approach  to  cointegration  is  tested  by  applying 


5  Results  of  lag  length  selection  are  available  from  authors  upon  request. 


Table  5 

Long-and-short  run  analysis. 


Dependent  variable : 

=  lnC02t 

Long  run  results 
Variable 

Coefficient 

Std.  Error 

T-Statistic 

Constant 

-50.2450 

4.7539 

-10.5692* 

lnYt 

10.0674 

0.8740 

11.5182* 

InY7 

-0.4942 

0.0410 

-12.0467* 

inCf 

0.8436 

0.0999 

8.43771* 

lnTRt 

0.0862 

0.0415 

2.07734** 

Short  run  results 
Constant 

0.0282 

0.0055 

5.0598* 

InY, 

4.9768 

1.6149 

3.0816* 

InY7 

-0.2392 

0.0807 

-2.9622* 

inCf 

0.4481 

0.0948 

4.7244* 

lnT/it 

-0.1090 

0.0453 

-2.4035** 

ECMt- 1 

-0.2994 

0.0948 

-3.1574* 

Diagnostic  tests 

Test  F-statistic 

Prob.  value 

x2  NORMAL 

0.8323 

0.6595 

X2  SERIAL 

1.9906 

0.1548 

x2arch 

0.0150 

0.9030 

X2  WHITE 

0.7055 

0.7110 

x2remsay 

2.3963 

0.1321 

*  Denote  the  significant  at  1%  level. 
**  Denote  the  significant  at  5%  level. 


Johansen  multivariate  cointegration  approach.  The  results 
reported  in  Table  4  confirm  that  there  are  two  cointegrating 
vectors  at  1%  level  of  significance  indicated  both  by  Trace  statistic 
and  Maximum  Eigen  values.  This  corroborates  the  long  run 
relationship  between  coal  consumption,  economic  growth,  trade 
openness  and  C02  emissions  in  case  of  India  over  the  study  period 
of  1966-2011. 

After  investigating  the  robustness  of  long  run  relationship 
between  these  variables,  the  long  run  marginal  impacts  of  coal 
consumption,  economic  growth,  trade  openness  on  C02  emissions 
in  case  of  India  is  investigated  by  applying  OLS  approach.  The 
results  reported  in  Table  5  postulated  that  coal  consumption  is 
major  contributor  to  energy  pollutants  and  it  is  statistically 
significant  at  1%  level  of  significance.  A  1%  rise  in  C02  emissions 
is  associated  with  0.8436%  consumption  of  coal.  These  findings 
support  to  the  view  by  Wolde-Rufael  [59]  that  coal  consumption 
highly  contributes  to  energy  pollutants  as  compared  to  other 
indictors  of  energy.  Jayanthakumaran  et  al.  [47]  also  nominated 
that  India  meets  55%  of  energy  demand  through  coal  consump¬ 
tion  and  produces  5%  of  C02  emissions  of  the  world6 7.  Further¬ 
more,  Jayanthakumaran  et  al.  [47]  reported  that  in  total  energy 
emissions  86%  is  contributed  by  coal  consumption  in  case  of  India 
and  same  inference  drawn  by  Alam  et  al.  [60]-7  . 

The  impact  of  trade  openness  on  energy  emissions  is  positive 
and  significant  at  5%  level.  All  else  is  same,  a  1%  expansion  in 
trade  openness  leads  to  upsurge  in  energy  emissions  by  0.0862%. 
This  finding  is  consistent  with  Jayanthakumaran  et  al.  [47]  also 
reported  positive  effect  of  international  trade  on  C02  emissions 
but  it  is  statistically  insignificant.  In  regional  studies,  Halicioglu 
[46]  reported  that  international  trade  increases  energy  emissions 
by  raising  economic  growth  and  same  inference  is  drawn  by  Nasir 
and  Rehman  [36]  in  case  of  Turkey  and  Pakistan,  respectively.  On 
contrary,  Shahbaz  et  al.  [37]  noted  that  trade  openness  improves 
the  environmental  quality  by  enhancing  the  capacity  of  country 


6  Other  sources  of  energy  such  as  crude  oil,  natural  gas  and  petroleum 
products  are  used  to  meet  the  rest  45%  energy  demand. 

7  Alam  et  al.  [60]  used  total  primary  consumption  as  an  indicator  of  energy. 
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Table  6 

VECM  Granger  causality  analysis. 


Dependent  variables 

Direction  of  causality 

Short  run 

Long  run 

ECTt_ , 

Joint  (Long-and-short  runs)  causality 

lnC02r 

lnYt 

InYf 

lnCt 

lnTm 

lnC02r,£Crt_ 

i  lnYt, ECT,_ 

l  lnYj'.ECT,., 

lnCt.ECTV, 

i  \nTRt,ECTt_ ! 

lnC02t 

_ 

3.6272** 

3.3713** 

12.8363* 

5.1278** 

-0.4550* 

- 

7.5935* 

6.9193* 

9.3793* 

7.8799* 

[0.0384] 

[0.0473] 

[0.0001] 

[0.0119] 

[-3.9339] 

[0.0006] 

[0.0040] 

[0.0001] 

[0.0005] 

<  m m Equation  id  =  ”eqn87” 

'/  4.6517** 

- 

8.1920** 

1.5550 

6.1098* 

-0.1632*** 

3.4984** 

- 

5.8590* 

1.2920 

4.5041** 

> 

[0.0174] 

[0.0000] 

[0.2227] 

[0.0060] 

[-1.8427] 

[0.0274] 

[0.0019] 

[0.2952] 

[0.0101] 

lnYj* 

4.2919* 

8.4460* 

- 

0.9678 

5.8392* 

-0.1535*** 

3.1922** 

6.0750* 

- 

1.1602 

4.3125** 

[0.432] 

[0.0001] 

[0.3878] 

[0.0072] 

[-1.7716] 

[0.0377] 

[0.0045] 

[0.3412] 

[0.01211 

InCf 

12.5233* 

0.5701 

0.5177 

- 

2.7016* 

-0.7704* 

10.0023* 

6.5615* 

6.6326* 

- 

7.2152* 

[0.0001] 

[0.5715] 

[0.6013] 

[0.0834] 

[-4.2526] 

[0.0001) 

[0.0015] 

[0.0012] 

[0.0009] 

In  Tri 

2.4178*** 

4.6856** 

4.5034** 

1.0847 

- 

- 

- 

- 

- 

- 

- 

[0.1078] 

[0.0169] 

[0.0195] 

[0.3509) 

*  Show  significant  at  1%  level. 

**  Show  significant  at  5%  level. 

***  Show  significant  at  10%, .level  in  parenthesis  we  reported  p-values. 


to  implement  advanced  technology  to  increase  domestic 
production. 

The  impact  of  linear  and  nonlinear  terms  of  real  GDP  per  capita 
is  positive  and  negative.  The  coefficients  of  linear  and  nonlinear 
terms  are  10.0674  and  -0.4942  and  it  is  statistically  significant  at 
1%  level  of  significance,  respectively.  The  significance  of  both 
linear  and  non-linear  terms  of  real  GDP  per  capita  provides  the 
empirical  evidence  of  inverted-U  shaped  relationship  between 
economic  growth  and  C02  emissions  “so  called  environmental 
Kuznets  curve”  (EKC).  The  empirical  exercise  revealed  that  a 
10.0674%  energy  emissions  are  linked  with  1%  increase  in  GDP 
per  capita  and  inverse  effect  of  squared  term  of  real  GDP  per 
capita  shows  the  delinking  point  of  C02  emissions  i.e.,  -0.4942 
and  real  GDP  per  capita,  once  an  economy  achieves  higher  level  of 
income.  This  empirical  evidence  provides  the  support  for  EKC 
revealing  that  C02  emissions  increase  at  the  initial  stage  of 
economic  growth  and  decline  after  a  threshold  point  i.e.,  Ind. 
Rupee  28,131  (US$531).  These  results  are  with  existing  literature 
such  as  He  [61],  Song  et  al.  [6],  Jalil  and  Mahmud  [20],  Fodha  and 
Zaghdoud  [22]  and  Jayanthakumaran  et  al.  [47]  for  China  and 
India,  Halicioglu  [46]  for  Turkey,  Lean  and  Smyth  [28]  for  ASEAN 
countries,  Nasir  and  Rehman  [36]  and  Shahbaz  et  al.  [37]  for 
Pakistan.  Saboori  et  al.  [65]  for  Malaysia  and  Shahbaz  et  al.  [68] 
for  Romania. 

The  short  run  dynamics  are  also  reported  in  lower  segment  of 
Table  6  and  the  results  indicated  that  coal  consumption  is  also 
main  contributor  to  C02  emissions  and  it  is  statistically  significant 
at  1%  level.  Interestingly,  international  trade  reduces  energy 
emissions  significantly  at  5%  significance  level.  The  coefficients 
of  linear  term  of  real  GDP  per  capita  and  nonlinear  i.e.,  squared 
term  of  real  GDP  per  capita  are  positive  and  negative,  respec¬ 
tively.  These  estimates  are  significant  at  1%  significance  level.  This 
further  confirms  the  existence  of  EKC  in  case  of  India.  The  short 
run  estimates  are  less  than  long  run  estimates  which  indicates  the 
reliability  and  stability  of  estimated  results. 

The  estimated  coefficient  of  the  lagged  ECM  term  is  -  0.2994  and 
significant  at  the  1%  level.  This  establishes  long  run  relationship 
among  the  running  variables.  This  suggests  that  changes  in  C02 
emissions  from  the  short  run  to  the  long  run  are  corrected  by  29.94% 
each  year.  The  significance  of  lagged  error  correction  term  further 
confirms  our  established  long  run  relationship  between  coal  con¬ 
sumption,  economic  growth,  trade  openness  and  C02  emissions. 

The  diagnostics  tests  such  as  LM  test  for  serial  correlation, 
normality  of  residual  term  and  white  heteroscedisticity  test  for 
short  run  model  are  reported  in  Table  5.  The  findings  suggest  that 


|  CUSUM  - 5%  Significance 


Fig.  3.  Plot  of  cumulative  sum  of  recursive  residuals.  The  straight  lines 
represent  critical  bounds  at  5%  significance  level. 

the  short-run  model  passes  all  diagnostic  tests  successfully.  The 
evidence  indicates  that  error  term  is  normally  distributed  and 
there  is  an  absence  of  serial  correlation.  There  is  no  evidence  of 
autoregressive  conditional  heteroscedisticity  and  same  inference 
for  white  heteroscedisticity.  The  short  run  model  has  passed 
Remsay  reset  test  which  confirms  that  the  functional  form  of 
short  run  model  is  correctly  specified. 

It  is  suggested  by  Browns  et  al.  [53]  that  stability  of  long  run 
and  short  run  parameters  can  be  tested  by  applying  CUSUM  and 
CUSUMsq  tests.  In  doing  so,  the  cumulative  sum  (CUSUM)  and 
cumulative  sum  of  squares  (CUSUMsq)  tests  have  been  employed 
for  parameter  stability. 

The  results  are  shown  in  Figs.  3  and  4.  The  graph  should  lie 
within  the  appropriate  bounds  if  the  parameters  are 
stable  (Pesaran  et  al.  [51]).  The  graph  of  CUSUM  at  the  5%  level 
and  we  can  claim  the  stability  of  long-and-short  runs  parameters. 


4.1.  The  VECM  Granger  causality  analysis 

The  presence  of  cointegration  among  the  variables  implies  that 
causality  relation  must  be  existed  at  least  from  one  side.  The 
directional  relationship  between  coal  consumption,  economic 
growth,  international  trade  and  C02  emissions  will  provide  help 
in  articulating  comprehensive  policy  to  economic  growth  by 
controlling  environment  from  degradation  and  utilize  energy 
efficient  technologies  importing  from  advanced  countries.  We 
applied  Granger  causality  test  within  the  VECM  framework  to 
detect  the  causality  between  the  variables.  Table  6  reports  the 
results  of  the  VECM  Granger  causality  analysis.  The  long  run 
causality  is  captured  by  a  significant  t-test  on  a  negative  coeffi¬ 
cient  of  the  lagged  error-correction  term  ECMt_r-  The  jointly 
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Fig.  4.  Plot  of  cumulative  sum  of  squares  of  recursive  residuals.  The  straight 
lines  represent  critical  bounds  at  5%  significance  level. 

significant  LR  test  on  the  lagged  explanatory  variables  shows 
short-run  causality. 

In  long  run,  our  empirical  evidence  reports  bidirectional 
causality  relationship  between  income  and  C02  emissions.  The 
feedback  hypothesis  exists  between  coal  consumption  and  C02 
emissions.  Coal  consumption  and  economic  growth  Granger  cause 
each  other.  This  finding  is  contradictory  with  Li  and  Li  [62]  who 
reported  that  coal  consumption  Granger  causes  economic  growth 
while  consistent  with  Paul  and  Bhattacharya  [63]  showed  feed¬ 
back  between  energy  consumption  and  economic  growth.  This 
suggests  that  government  should  encourage  energy  exploration 
polices  which  in  resulting  spur  economic  growth  and  hence 
demand  for  energy.  There  is  bidirectional  causal  relationship 
between  economic  growth  and  C02  emissions.  These  findings 
confirm  our  long  run  as  well  short  run  results  that  income  has 
positive  and  negative  effects  on  C02  emissions.  This  further 
validates  the  existence  of  environmental  Kuznets  curve  (EKC). 
On  contrary,  Ghosh  [26]  and  Alam  et  al.  [60]  presented  no  causal 
relationship  between  income  and  C02  emissions. 

This  shows  the  significance  of  our  empirical  evidence  as  we 
applied  more  specified  economic  model.  Trade  openness  Granger 
causes  economic  growth,  coal  consumption  and  C02  emissions. 
This  implies  that  trade  openness  stimulates  economic  activity  and 
adds  in  energy  demand  (coal  consumption)  which  emits  C02 
emissions.  In  short  run,  bidirectional  causality  is  found  between 
economic  growth  and  C02  emissions,  coal  consumption  and  C02 
emissions,  trade  openness  and  C02  emissions.  The  findings 
indicated  that  the  long  run  income  elasticity  for  C02  emissions 
is  less  than  the  short  run  elasticity  is  further  evidence  for  an  EKC 
(see  Narayan  and  Narayan  [30]  for  more).  The  feedback  hypoth¬ 
esis  also  exists  between  trade  and  economic  growth.  The  unidir¬ 
ectional  causal  relationship  is  found  running  from  trade  openness 
to  coal  consumption  and  no  causality  exists  between  coal  con¬ 
sumption  and  economic  growth.  The  results  of  joint  causality 
confirm  the  long-and-short  runs  findings. 


5.  Conclusion  and  policy  implications 

This  paper  investigated  the  dynamic  linkages  between  coal 
consumption,  economic  growth,  trade  openness  and  C02  emis¬ 
sions  in  case  of  India  over  the  period  of  1966-2011.  We  have 
applied  Ng-Perron  [57]  unit  root  test  to  test  stationarity  of  the 
variables  and  Narayan  and  Popp  [49]  is  accompanied  to  examine 
the  robustness  of  empirical  evidence.  Further,  the  ARDL  bounds 
testing  approach  to  cointegration  has  employed  to  investigated 
the  long  run  relationship  between  the  variables  and  Johansen 
multivariate  cointegration  is  applied  to  test  the  strength  of  long 
run  results.  The  present  study  has  important  policy  implications 
for  Indian  economy.  The  empirical  evidence  of  strong  relationship 
between  coal  consumption  and  economic  growth  indicates  that 


sustained  level  of  economic  growth  can  be  achieved  by  exploring 
the  new  sources  of  energy  supply  for  long  span  of  time. 

Our  results  confirmed  that  all  the  variables  are  cointegrated 
for  long  run  relationship  and  long  run  results  are  robust.  Coal 
consumption  is  main  contributor  to  C02  emissions.  A  rise  in 
economic  growth  raises  demand  for  coal  that  emits  energy 
pollutants.  Trade  openness  is  also  positively  linked  with  C02 
emissions  but  in  short  run  former  declines  later.  Results  also 
validated  the  existence  of  environmental  Kuznets  curve  (EKC)  in 
long  run  as  well  as  in  short  run  in  case  of  India.  The  causality 
analysis  explored  the  bidirectional  causal  relationship  between 
economic  growth  and  C02  emissions  and  feedback  hypothesis 
exists  between  coal  consumption  and  C02  emissions.  Coal  con¬ 
sumption  and  economic  growth  Granger  cause  each  other.  Trade 
openness  Granger  causes  economic  growth,  coal  consumption 
and  C02  emissions.  Economic  growth  and  C02  emissions  Granger 
cause  trade  openness  in  short  run. 

Hence,  our  results  shows  that  dependence  on  the  coal  con¬ 
sumption  as  a  factor  that  promotes  economic  growth  might  be 
feasible  option  for  short  run  but  in  the  long  run  we  have  to  look 
for  alternative  of  this.  This  is  particularly  because  coal  consump¬ 
tion  increases  C02  emissions  and  which  indirectly  has  negative 
impact  on  the  economic  growth.  The  other  interesting  issues  are 
the  Granger-causal  relationship  between  C02  emissions,  trade 
openness,  and  coal  consumption.  Since  trade  openness  is  found  to 
Granger-cause  C02  emissions,  we  can  suggest  Indian  policy¬ 
makers  on  a  few  important  points:  (1)  allowance  of  foreign 
investors  to  invest  in  the  India  particularly  in  know-how  and 
advanced  technologies  (2)  free  entry  of  FDI  because  when  FDI 
flows  into  a  host  country,  FDI  may  transfer  new  ideas,  technol¬ 
ogies  and  managerial  skills  to  domestic  firms  and  which  in  turn 
expected  to  reduce  C02  emissions  (3)  promotion  of  such  activities 
through  which  foreign  firms  may  provide  continuous  technologi¬ 
cal  assistance,  such  as  introduction  of  new  technologies  and 
training  of  workers,  to  the  domestic  firms  which  will  lead  to 
decrease  in  C02  emission  in  one  hand  and  increase  economic 
growth  on  the  other.  These  processes  are  also  expected  to  reduce 
dependence  on  the  consumption  of  coal  without  harming  the 
interest  of  higher  economic  growth. 
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